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Abstract. We formulate deformation of relativistic stars 
due to the magnetic stress, considering the magnetic fields 
to be perturbations from spherical stars. The ellipticity 
for the dipole magnetic field is calculated for some stel- 
lar models. We have found that the ellipticity becomes 
large with increase of a relativistic factor for the models 
with the same energy ratio of the magnetic energy to the 
gravitational energy. 
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1. Introduction 

There is a growing interest in new classes of objects: soft- 
gamma repeaters (SGRs) and anomalous X-ray pulsars 
(AXPs). Until now, four identified SGRs are reported ob- 
servationally. Their associations with supernova remnants 
strongly suggest that the SGRs are young neutro n star s 
(see e.g. Kulkarni & Frail |1993[ Murakami et al. [19941) . 



Furthermore, recent measurements (see e.g. Kouveliotou 
et al. 1998 , 1999) of the period and period derivative 



yield evidence for these pulsars to be ultramagnetized 
neutron stars with field strength (~ lO^^G) in excess of 
Bci ^ lO^^G, i.e., magnetars (Duncan & Thompson 1993; 



Thompson & Duncan 1993, 1995, 1996). Some class of 
X-ray pulsars also suggests the magnetic fields with 10^''- 



lO^^G (see e.g. Mereghetti & Stella |1995D . Such magnetic 
fields are much stronger than that of known pulsars (10^- 
lO^^G; see e.g. Taylor et al. |1993D until then. Though the 
relation between the SGRs and the AXPs is not yet clear, 
there exist neutron stars with very strong magnetic fields. 
In these ultramagnetized stars, the magnetic influence be- 
comes important as well as the relativistic effects. If we as- 
sume that a long-lived electric current flows in highly con- 
ductive neutron-star matter, the magnetic pressure corre- 
sponding to the Lorentz force conies into play. Hence, it 
induces deformation of stars. In this paper, we study such 



deformation from spherical stars within a general relativis- 
tic framework. 

The quadrupole deformation of magnetized Newtonian 
stars was discu ssed by Chandrasekhar & Fermi ( 1953 ) and 
Ferraro ( 1954 ), in which incompressible fluid body with a 
dipole magnetic field is assumed. This deformation has 
been discussed also in rel ation to the gravitational radia- 
tion (Gal'tsov et al. |l984t Gal'tsov & Tsetkov |l984|) The 
general relativistic approach by B ocquet et al. ( |1995| ) and 
Bonazzola & Gourgoulhon (1996) has appeared recently. 
However, their approach is fully numerical. In this paper, 
we develop almost an analytical treatment by assuming 
weak magnetic fields compared with gravity. This assump- 
tion is valid even in the magnetars. Our formulation is re- 
garded as a general relativistic version of Chandrasekhar 
& Fermi (1953) and Ferraro (1954). In our method, we 
can easily include realistic equations of state (EOS) and 
construct relativistic magnetized stars. Furthermore, this 
method gives simple calculations of ellipticity of deformed 
stars, and so on. 

Since the observed ultramagnetized neutron stars have 
long periods (T ~ several sec), we may neglect the rota- 
tion of the magnetized stars, that is, we discuss static 
cases. We take non-rotating, spherical relativistic stars as 
backgrounds, and consider the magnetic fields as the per- 
turbation. In particular, we consider only axisymmetric, 
poloidal magnetic fields produced by long-lived (toroidal) 
electric currents, because toroidal magnetic fields would 



break the symmetric property (see also Bocquet et al. 1995| 
and reference therein). Furthermore, we assume a per- 
fectly conducting interior. Since we now consider non- 
rotating configurations, this implies that the electric field 
inside the stars must be zero. Hence, there is no electric 
charge inside the stars. From this, we can write the 4- 



current as J^ = (0, 0, 0, J^) (Bocquet et al. 1995). Further- 



more, the surface charge should be absent, since the total 
charge should vanish in astrophysical situations. Other- 
wise, the electromagnetic field itself would have the an- 
gular momentum due to the non-vanishing electric field 
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produced by th e charge (Feynman et al. |1964| ; Ma |1986 
de Castro 1991). This is not the purely static case. 
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The current distribution is introduced as the first-order 
quantity with respect to the perturbation. The corre- 
sponding magnetic field is solved by the Maxwell equation. 
We shall investigate deformation of stars due to the resul- 
tant magnetic stress, which arises as the second-order cor- 
rection to the background field. This perturbation method 
is very sim ilar t o that of slowly rotating stars developed 
by Hartlc ( |l967|) , in which the rotation is regarded as a 
small parameter. Our formalism can be applied to any 
configurations of the magnetic fields. However, we restrict 
ourselves to dipole magnetic fields because the dipole fields 
are important in many astrophysical situations. 

The plan of this paper is as follows. In Sect. 2, the mag- 
netic fields are investigated in the background space-time. 
The effect arising from the magnetic stress on equilibrium 
of stars is considered in Sect. 3. The solution correspond- 
ing to the quadrupole deformation of the stars will be 
given. To evaluate the deformation quantitatively, the el- 
lipticity is discussed in Sect. 4. Finally, Sect. 5 is devoted 
to the discussion. Throughout this paper, we use the units 
in which c = G = 1 , and the Gaussian system of units for 
electromagnetic fields. 



2. Magnetic fields in spherically symmetric space- 
time 

We now consider an axisymmetric, poloidal magnetic field 
created by a 4-current J^ = (0, 0, 0, J^) in a non-rotating, 
spherical star. We suppose that the magnetic field is weak, 
i.e., Bf^ ~ 0{e). The line element of the spherically sym- 
metric space-time is given by 



ds2 ^ -e'^^'-Uf + e^f'^'dr^ + r^dd^ + r^ sin^ Od^j)^, (1) 



where v and A are functions of the radial coordinate r only. 
Since we do not consider electric fields, we can assume 
that the electromagnetic 4-potential Afj_ has only the (p- 
component, i.e., A^ — {0,0,0, A^) (see also Bocquet et 
al. 1995). In this case, the Maxwell equation is reduced to 



,-A 



Qj.2 



1 



+ -{v'~ A') , 



.dA^ 
dr 



I d^As 1 ^„dA 

-rr cot 



r2 96*2 



de 



-47rJrf, 



(2) 



where the prime denotes the differentiation with respect 
to r. 

We expand the potential A^ and the current J^ as 
follows (Regge & Wheeler 1957): 



A 



2_^ai{r)sm0 — - 

1=1 

J4> = 2^3i{r)s\ne 



1=1 



de 



(3) 
(4) 



where Pi is the Legendre's polynomial of degree I. Substi- 
tuting these forms into Eq. (2), we have 



-x(Pai 



dr^ 



— '+^(.'-A')e-$ 



l{l + l) 



dr 



-ai 



-Anji. (5) 



For a given current ji, we can obtain the potential ai and, 
therefore, the magnetic field. From now on, we only con- 
sider a dipole magnetic field, i.e., / = 1. The potentia l out- 
side the s tar is easily solved (Ginzbur g fc O zernoi 1965 ; 
Petterson 1974 ; Wasserman & Shapiro |1983|) in the form 



fli 



3/i_^2 



8AP 



In 1 



2Af 



2Af 
r 



2M2 



(6) 



where /x is a constant corresponding to the magnetic dipole 
moment with respect to an observer at infinity, and M is 
the total mass of the background star. In order to de- 
scribe the magnetic field inside the star, we require the 
current distribution ji . The current ji is not arbitrary but 
subject to an integrabilit y cond ition (Ferraro 1954; Ghan- 
drasekhar & Prendergast 1956; Bonazzola et al. 1993). As 
will be shown in Eq. (Eq), this current is given, up to the 
first order in e, by 



jiir) = cor^ iPo{r) + Po{r)) 



(7) 



where cq is an arbitrary constant, and po and po denote 
the density and pressure, respectively, of the background 
star. By requiring that oi behaves as a regular function 
at the center of the star, we now obtain the potential ai 
in the vicinity of the center: 



ai 



aor 



0\ 



(8) 



where ao is a constant, which is fixed by the boundary 
condition at the surface. In this way, we can construct the 
magnetic field in the whole space-time. 

Fig. |l| displays the tetrad component of the magnetic 
field, 



H 1 « A 2cos^ 

'■ ^ 2~- — 7)'^^^"^ ~ 2 '^1 



(9) 



on the symmetry axis (i.e., 9 — 0), and Fig. g displays the 
tetrad component 



Ba 



rsin^ 



drA^ 



(10) 



on the equatorial plane (i.e., 9 — 7r/2) with respect to 
the radial coordinate r. We have normalized them by the 
typical magnetic field strength ii/ B? , where B is the ra- 
dius of the star. The sohd hues denote a relativistic case, 
whereas the dashed lines correspond to a Newtonian case. 
In these calculations, we have used the polytropic EOS: 
Po = kPq (7 = 2). From these figures, we see that the in- 
tensity of the magnetic field increases as r becomes closer 
to the center. Furthermore, these two figures show that 
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Fig. 2. The tetrad component of the magnetic fields, Bg on the 
Fig. 1. The tetrad component of the magnetic fields, Bf on equatorial plane (6» = 7r/2), plotted against r/R. The solid Une 
the symmetry axis (6^0), plotted against r/R. The solid hne denotes a relativistic case {M/R = 0.2), while the dashed hne 
denotes a relativistic case {M/R = 0.2), while the dashed line corresponds to a Newtonian case {M/R = 0.01). The magnetic 
corresponds to a Newtonian case {M/R = 0.01). The magnetic ggi^s are normahzed by the typical field strength ^^/R\ 
fields are normalized by the typical field strength n/R . 



despite the same magnetic inom.ent with respect to an 
observer at infinity, the central magnetic field of the rel- 
ativistic star is stronger than that of the Newtonian star 
by about 50% of the Newtonian case. Therefore, it follows 
that the relativistic effect strengthens the internal mag- 
netic fields. 



The total energy-momentum tensor is the sum of the 
perfect-fluid part Tj^,,** and the electromagnetic part 



J- (cm) , 






where 



(12) 



3. Equilibrium configurations of magnetized stars 

Next, we consider deformation of magnetized stars due 
to the magnetic stress, which is regarded as the second- 
order effect. We formulate the deformation of the star and 



space-time, following Hartle ( |1967D 



r„ 



{p + p)u^'u^+p5^'^, 



T^. - i {f'-'F'^^ - \f.^f"'k) ■ 



(13) 
(14) 



In Eq. (14), F^i, is the Faraday tensor. The pressure p and 
the energy density p can also be expanded in multipoles 
as 



3.1. Equations of equilibrium 

The metric can be expanded in multipoles around the 
spherically symmetric space-time. In particular, when we 
deal only with a dipole field, i.e., Z = 1 in Eqs. (3) and 
(4), t he metric can be written i n the form (see also Hartle 



p{r, 0) = po + {Sp{i=o) + 5p(i=2)P2) , 

p(r, 0) ^ po + -y ((5p(,=o) + Sp{i=2)P2) , 
Po 



(15) 
(16) 



1967| ; Chandrasekhar & Miller |1974| ) 



where (^P(;=o) ^-nd (5p(;=2) depend on r only, and we have 
assumed a barotropic case. 

From the Einstein equation, we can obtain 



ds^ = -e''^''^[l + 2{ho{r) + h2ir)P2{cos0))]df 



=.A(r) 



(mo(r) +m2{r)P2{co8 0)) 



dr' 



+ r^[l + 2k2{r)P2{cos0)] {d0^ + r^ sin^ 0d(f>^) ,{U) 

where ho, ft.2, "^Oj "I2 and ^2 are the corrections of the 
second order in e. 



mo = 47rr2-^(5p(i=o) + x 



-,-A/„m2 



ia[r + ^ai 



Ii'q ~ Anre Sp^i^^) + -^' ^ "^0 + ^e mo 



e 
3^ 



-A 



{<r 



(17) 



(18) 
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ft, + Kt = "-2 77 H "i2 

\r 2 / r 






^ 4 , 




(19) 


h2 + -^2 - ?e-^ (al)' , 




(20) 



h' 



Ae^ 



7W2 



7r-(po+Po) + -5-7(l-e')-'^' 



Sr'^iy 



■e^al + 



3r3.'V' + ~'"^"^ 



/I / A 2 \ , ,,2 167r > 



(30) 



where V2 = /i2 + fc2- These equations govern the relativistic 
magnetized star. We have to solve four differential equa- 
tions (27)-(30) and one algebraic equation (20) for the 
unknown functions {mQ,m2,ho,h2,V2). Furthermore, we 
can derive 6p(^i=o) and Sp(i^2) by substituting the solution 
of Hq and /i2 into Eqs. ( [2q ) and (23). 
(21) In order to solve Eqs. (27) and (28) inside the star, it 

is also convenient to introduce a quantity 
Furthermore, from the conservation law of the total 



^2 + ^2 + 2'^'^' k'2 

= 47rre^(5p(i=2) + — e^m2 + -v'e^m2 



3a, 2 ., 1 > 

r r 6r 



e-\a[f + ^al 



energy-momentum tensor, we obtain 



-i.'(^ + i)%.., 



Sp{i=2) = - {PQ + Po) h2 ~ ^aiji, 

2 
- (Po+Po)/i2" ■^"i-?i- 



po+po 
From Eq. (p6|), we have 



(22) 2 



-ai = ci. 



(31) 



(32) 



(23) 



(24) 



The integrability condition for Eqs. (22) and (24) leads to 



r'^{po+Po) 



const. (= Co) . 



(25) 



3r^{po+Po) 
Moreover, Eqs. (27) and (28) are rewritten as 

m', = 47Tr^Pk^p„+p„)5P^ + y-^(^a[f + l_al (33) 

SPq = - Unpo + ^ j e^^TTiQ - Anre^ {po + Po) SPq 

1 / /n2 ,2 ;^ 2 2 jl 

dr 3r-^ 3 r-^ [po + po) 



a[. (34) 



In the next subsection, we solve these differential equa- 
Thisis consistent with Eq. (5.29) of Bonazzola et tions for the metric functions, 
al. (1993) up to the first order in e. Using this current 
distribution, Eq. (22) can also be integrated as 3^2. The exterior solution and boundary condition 

Sp^i^o) =-{po+ Po) hn + ^aiji + ci (po + Po) , (26) First, we consider the solution outside of the star, in which 

3'' Pq = Pq = and ji — 0. 



where ci is a constant of integration. 

Consequently, we have two set of differential equations. 



TOn 



47rr2^(po+ Po)(/io-ci) 



mo 



Pq 

-,-A/„m2 



The solution of too and /lo is given by 

2A/ 



%M^r 



3^ (J2. ,f„ ,,r2\ 



871 p^, 



+ 3e-^Kr + ^«f + y-aUt, 
h'fy = [ — H ) e^mo - Anre^ (po + Po) (ho - ci) 



,2 
1 
3^ 



+ - K) - 3,3 



2 2 A„ 2 Stt , 



e ai + -—e aiji, 
6r 



(27) 



(28) 



8M5 



r^ - Mr -Ar In 1 



In 1 



and 



+ -^r^ (r ~ 2Af ) 
32Af 6 ^ ' 

3/i2 4r - M 
~ 8M3 r (r - 2A-/) 

3^2 (r-Af)(r-3A/) 
"*" 8Af 5 r-2A/ 



2Af 



+ C2, (35) 



~z^'/i2 + -e 



2 -A /I i^' 



+ ^)Wif + ^^i< (29) 



3p2 
32A/6 ' 



In 1 



2Af 



ln( 1- 

2 



2Af 
r 

C2 



r - 2M 



C3, (36) 
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where C2 and C3 are constants of integration. At large r, 
mo and ho behave as 



mo 
ho 



C3 



^4 

3/i^ 

8M4 



O 



1 



£2 
r 



O 



1 



Since ho must vanish at infinity, we obtain 



C3 



8M4' 



The constant ci corresponds to the mass shift, which is 
fixed by matching with the interior solution at the surface. 
On the other hand, the differential equations for v^ 
and /i2 is rather complicated, but analytically solved. The 
solution of V2 is 



V2{z) 



IK 



V^2zri 



Q\{z)^ 



3/i2 



4A/4Vz2 - 1 



P2'(z) + t;2pW,(40) 



where z is defined as 2: = r jM — 1, iiT is a constant of inte- 
gration, P\ and Q\ are the associated Legendre functions, 
and vip is 



W2p 



V 3m^ 7z2-4 



4M4 8M4 z2 - 1 

3/i2 Z (llz2 -7) ;2 _ 1 

^r m 



16M4 
,2 



16M4 ^ '' V 2 + 1/ 



Furthermore, /12 is given by 
3Ai2 



h^{z) = A'Q^(z) 



8Af' 



P2'W+/i2pW 



(41) 



(42) 



where P| and Q\ are the associated Legendre functions, 
and /i2p is written as 



hip — 



16M4 

3m' 
32M4 

16Af4 



6z^ 



3z - 6 



4z^ + 2z' 



1 



32^-8 



{^ - 1) 



z-3 

In 



z- 



In- 



z- 1 

z + 1 



z + 1 



(43) 



In Eqs. (^) and ([4^), we have used the boundary condi- 
tion at infinity. The remaining constant K will be fixed 
by the boundary condition at the surface. Thus we have 
obtained the analytical solution outside the star. 

We turn our attention to the interior solution. For a 
given EOS, we can obtain the solution numerically. For 
the actual numerical work, we investigate the behavior of 
the metric functions in the vicinity of the center. 



First, we consider the metric functions mo and bPo- 
The solution in which both mo and 8P0 vanish at the 
center (see also Chandrasekhar & Miller 1974) is given by 

(44) 

(45) 

The regular solution at 



(46) 

where constants /3i and /?2 are not independent by the 
regularity condition at the center. 

Finally, we can obtain the metric functio ns by impos- 
ing the junction conditions (O'Brien & Synge 1952 ) at the 
surface: 



(37) 


2 2 3 


(38) 


(5^0 - -0 ("0^ ~ coao) r' H 




where ao is defined in Eq. (8). 




Next, we consider v^ and hi 


(39) 


the center is 


vi ~ /?ir4 -)-•••, 


;h is 


hi ~ Pir' + • • • , 



9i^A+n = di^A-R {n,iy = t,r,( 



(47) 
(48) 



where g^^ denotes the metric components. From these 
conditions, the integration constants Ci, K, Pi and j3i are 
fixed. 

4. Ellipticity of magnetized stars 

We consider the magnetic field on the stellar shape of the 
equilibrium. The additional Lorentz force J x B mainly 
acts on it in the perpendicular direction to the symmetry 
axis {6 = 7r/2), that is, fiattens the star. The flattening ef- 
fect is also recognized by considering the (r, r) component 
of the magnetic stress tensor, 

T,.n/, = ^{BeB<'-B^B^). (49) 

Along the symmetry axis, Bg must vanish owing to the 
axisymmetry. Hence, the stress Tf^^^-f is negative on this 
axis. On the other hand, on the equatorial plane, since 
Br is zero at any r, the stress T'(cm)'^ has the opposite 
sign. This indicates that the spherical star is shrunk in 
the parallel direction to the symmetry axis {0 — 0) and 
expanded in the perpendicular direction {9 = 7r/2) by the 
magnetic effect. Thus we can see the flattening effect. 

Next, in order to evaluate the deformation quantita- 
tively, let us introduce the ellipticity, which is defined as 



ellipticity 



_ (equatorial radius) — (polar radius) 



(50) 



(mean radius) 

where these radii denote the circumferential radii under 
general relativistic situations. From this definition, ellip- 
ticity is given by (Chandrasekhar & Miller 1974) 

u- +■ ■-^ 2co 3hi 3 
ellipticity ~ — -ai -\ ; — -/c2. (51) 
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The first term of the right-hand side in Eq. (51) corre- 
sponds to, in a sense, the effect of the 'Lorentz force', 
the second term represents the effect of the perturbation 
of the 'gravitational potential' induced by the magnetic 
effect, and the third term is a 'purely relativistic term' 
which arises from a definition of the radius, that is, the cir- 
cumferential radius. These contributions to the ellipticity 
are shown in Fig. || as a function of the relativistic factor 
M/R. We have normalized them by jjl^R^ /P [I is the mo- 
ment of inertia, which is well defined also in the relativistic 
case), and we have used a polytropic EOS (7 = 2). From 
this figure, we can see that the term due to the 'gravita- 
tional potential' does not change significantly, while the 
term concerning the 'Lorentz force' increases with the rel- 
ativistic factor M/R as known from Figs. |l| and |[ 



Fig. 4. The variation of the ellipticity with respect to M/R. 
The solid line denotes the case of the polytropic EOS with 
7 = 2, and the dashed line corresponds to the case of 7 = 5/3. 



culations give the same results (see Fig. 2 of their paper), 
that is, (3 takes the values very close to 1. 



Fig. 3. The three contributions to ellipticity: (a) the effect of 
the 'Lorentz force', (b) the effect of the perturbation of the 
'gravitational potential', and (c) the 'purely relativistic effect' 
(see text). These are plotted against M/R. We have used the 
polytropic EOS (7 = 2). 



Fig. H displays the (total) ellipticity for different poly- 
tropic models (7 = 5/3 and 7 = 2). From this figure, 
we find that the ellipticity becomes large as the relativis- 
tic factor M/R increases, in each case of 7 = 5/3 and 
7 = 2. The common feature of the monotonic increase 
shows the effect of the 'Lorentz force term' to be effec- 
tive. An important thing is that the relativistic calcula- 
tion leads to much larger ellipticity for a fixed value of 
/i^i?^//^. Finally, we give a comment concerning the pre- 
vious general-relativistic studies. The quantity plotted in 
Fig. y is exactly the magnetic distortion factor f3 intro- 
duced by Bonazzola & Gourgoulhon ( 1996 ). It is noticed 
that although the EOS which we have used is different 
from that of Bonazzola & Gourgoulhon (1996), both cal- 



5. Discussion 

Recent observations suggest that some class of neutron 
stars has strong magnetic fields. These may promote a 
new branch with magnetized stars. As a simple approach 
to the models, we have formulated the structure of the 
magnetized stars within a general relativistic framework, 
considering the perturbation from a spherical star. In par- 
ticular, a dipole magnetic field has been dealt with. We 
have showed the current distribution which yields equi- 
librium configurations up to the second order in e. Fur- 
thermore, the ellipticity of the star has been estimated as 
a simple example. We have found that the ellipticity be- 
comes large as the relativistic factor M/R increases, for 
the same energy ratio of the magnetic energy to the grav- 
itational energy. Our analytical approach have made the 
calculations much simpler than that of the previous work. 
This method can be extended to more general cases of re- 
alistic EOS and general current distribution, in which the 
current exists in some domain of the star. Therefore, this 
can be applied to wider range of astrophysical situations. 
Another extension of this work is to incorporate ro- 
tation of stars. The stationary configurations, in which 
the rotation axis is aligned with the magnetic axis, make 
the calculations complex because of appearance of non- 
vanishing electric fields. However, this can be managed. 
Since the rotational effect deforms the star as well as the 
magnetic effect, we are also interested in seeing which of 
them to be effective. These will be the subject of further 
investigation. 
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